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Ensemble averaged coherent state path integral for disordered bosons 
with a repulsive interaction and derivation of a nonlinear sigma model 

Bernhard MieckQ 
Abstract 

A coherent state path integral is considered for bosons with an ensemble average of a random potential and 
with an additional, repulsive interaction in the context of BEC under inclusion of specially prepared disorder. The 
essential normalization of the coherent state path integral, as a generating function of observables, is obtained from 
the non-equilibrium time contour for 'forward' and 'backward' propagation so that a time contour metric has to be 
taken into account in the ensemble average with the random potential. Therefore, the respective symmetries for the 
derivation of a nonlinear sigma model follow from the involved time contour metric which leads to a coset decomposition 
Sp(4) /U(2) (gi U(2) of the symplectic group Sp(4) with the subgroup U(2) for the unitary invariance of the density- 
related vacuum or ground state; the corresponding spontaneous symmetry breaking gives rise to anomalous- or 'Nambu'- 
doubled field degrees of freedom within self-energy matrices which are finally regarded by remaining coset matrices. 
The notion of a 'return probability', according to the original Anderson- localization', is thus naturally contained within 
coherent state path integrals of a non-equilibrium contour time for equivalent 'forward' and 'backward' propagation. 
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1 COHERENT STATE PATH INTEGRAL AND AVERAGING METHOD FOR DISORDER 



1 Coherent state path integral and averaging method for disorder 

1.1 Ensemble averages for static and dynamic disorder 

The original concept of 'Anderson localization' is combined to a 'return probability' of a wave packet, which starts out 
to propagate from an initial space point and which is scattered back by the impurities and the disordered potential, for 
a measure of its 'localization' [T]. This 'return probability' can therefore be specified by the propagation of ensemble 
averaged, retarded and advanced Green functions [2 . Our presented approach in terms of coherent state path integrals, 
adapted to a system of BEC in an external trap potential, naturally contains this concept by a time contour metric with 
'forward' and 'backward' propagation^- [T3] . The ensemble average of a random potential involves the time contour metric 
which thus determines the symmetries for the derivation of a nonlinear sigma model [T51 116] . The coset decomposition 
Sp(4) / U(2) (g) U(2) of an anomalous-doubled self-energy gives rise to a spontaneous symmetry breaking (SSB) with 
the subgroup U(2) for an invariant vacuum or ground state of density-related field degrees of freedom. The remaining 
coset matrices Sp(4) /U(2) finally comprise the anomalous- or 'Nambu'-doubled field degrees of freedom from the off- 
diagonal block parts of the total self-energy which is accomplished by a Hubbard-Stratonovich transformation (HST) of 
anomalous-doubled, dyadic products of boson fields with incorporation of 'hinge' fields from the SSB [T71 fT51 [T5] . 

We consider two models (jl.ip with Hamiltonians Hi{^ , -0, Vj) and Hijffl, Vu) in normal ordering of bosonic 
creation and annihilation operators and with static and dynamic disorder potentials Vi(x), Vu(x,t) in parallel. Both 
Hamiltonians (jl.ip contain the same one-particle part h(x) (|1.2I) with the kinetic energy and with an external trap- 
potential u(x), which is shifted by a reference energy /j, as a chemical potential. Apart from the same one-particle part 
h(x) (jl.2l) . we take the identical repulsive interaction with parameter Vq > for a quartic contact interaction of bosonic 
operators '^t & ipg ipg Furthermore, even- and complex-valued, spatially local source fields j^-sif), ji)>ip;x{t) are 

included in common in order to allow for a SSB with a coherent, macroscopic wavefunction and with anomalous (or 
'Nambu') paired bosons, respectively. These source fields can also be used for the determination of observables from 
differentiating various manners of the generating function. According to a presupposed, spatially spherical symmetry of 
u(x), we normalize the spatial summations Y]g . . . (| 1 . 3|) by the spherical system volume of '£)' dimensions and by the 
spatial unit cell (ax) d which yields a parameter M x = /(ax) d to be applied for various approximations, as e.g. in a 
saddle point computation. The time parameter t p is restricted to the range . . . + To with the discrete intervals At > 
so that one is limited by the order of maximum energy h£L max = K/ At within the propagation of the two Hamiltonians 
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The two Gaussian distributions for random potentials Vi(x), Vu(x,t) are determined by the second moments (jl.4|1.5|l for 
static and dynamic disorder, respectively, and vanishing mean values. Both distributions are delta-function like, concern- 
ing the spatially 'contact' Kronecker-delta and concerning the white-noise delta-function of time. Moreover, we emphasize 
the two different normalizations of second moments in (|1.4I1.5|) which are important in subsequent transformations and 
derivations for a proper, finite scaling of energy ranges within the nonlinear sigma models (cf. the second moments 
and their normalization in random matrix theories). Therefore, there occur two different disorder parameters Ri, Rn 
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of physical dimensions [energy x time] and [energy x (time) ] in the two moments (|1.4I1.5[) . We apply the normalized 
generating functions Z[3,Vi], Zu[3,Vu] (|1.6p with the time development operators (11.71) which are composed of the 
prevailing Hamiltonian T or 'II', each having an additional source matrix 3 for calculating observables by differentiating 
the coherent state path integral. The latter path integral (|1.6I1.7I) results into proper normalization of unity (|1.8[) for 
vanishing source matrix 3 = and identical values (|1.9p of symmetry breaking fields jib- t 3(t), j^-sit) on the two branches 
of propagating time development operators in forward and backward direction 
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In accordance to the forward '+' and backward '-' propagation of time development operators, we introduce a contour 
time '£p=±' and contour integral J c dt p . . . with a contour time metric , n p= ± = ±1' in order to regard the changing sign of 
phases in the exponent of the time evolution operators (|1. 11111. 14| [20]- [27]. In the remainder we will also briefly outline 
the various steps of transformations to a nonlinear sigma model with incorporation of the precise time steps and intervals 
within the coherent state path integral which is specified by the normal ordering of creation and annihilation operators. 
These precise time steps, with additional time shift 'a£ p ' in the complex conjugated coherent state fields, are usually 
omitted for brevity in literature, but are ubiquitous in many-particle physics; this problem of proper time steps does 
not prevent transformations and derivations to sigma models with SSB and a coset decomposition and has necessarily 
to be regarded as soon as quantum mechanical field integration variables at neighbouring, but still different times 
and Hp + At p \ are considered within normal ordered coherent state path integrals. We briefly hint at the problem of the 
appropriate, precise time steps by defining a slightly modified time contour (|1.13I1.14[) for simplified representation of the 
exact time steps (cf. also additional boundary conditions of coherent state fields in (j2.22M2.25)) ) 



dt 



p ■ ■ 



+T 



dt A 



dt 



p ■ ■ 



E 



P =± jo 



-T 



dt T . 



r+T y.+T 

dt- ■■■= I dt+ ■■■ dt- 

+T JQ JO 

... ; n p = = +1 ; ?7_ = -1 j ; 



dt,, 



dt,, 



+T 



dh 



-At 



+ (T +At) 



dt- 



+T M 

dt + ...+ dt 

-At J+(T +At) 

_ r+* P 

22 I dt p Ti p ... ; 0+ = -At ; T+ = T ; 0_ = ; T_ = T + At . 
P =± J?>p 



(1.11) 
(1.12) 

(1.13) 
(1.14) 



Since coherent state path integrals allow for the exact time sequence of coherent state fields with proper, additional 
time shifts 'At p ' of the corresponding complex-conjugated fields, one can also investigate other kinds of coherent state 
path integrals as those of Eqs. (11.4111.10]) [281 [29l [30] . Apart from the presented problem, we have also applied coherent 
state path integrals to a trace representation of delta functions with maximal commuting sets of symmetry operators 
(|1.15ll.l6p . As these symmetry operators (| 1 . 1 7fl 1 . 2U]l are given in terms of normal ordered creation and annihilation 
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2 ENSEMBLE AVERAGES IN MODEL I AND II FOR A NORMAL-ORDERED HAMILTONIAN 



operators, one achieves a similar coherent state path integral on time contours as (11.4111. 1011. 11111. 14")) 
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Coherent state path integral representation of the trace 
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:= maximal commuting set of second quantized operators according to symmetries; (1-16) 
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This is accomplished by the application of the Dirac identity (|1.21[) to the principal value '*P' and delta function of the 
symmetry operators (| 1 . 22ti 1 . 25]) so that propagation with exponentials is also implied on two branches of a 'disconnected 
time contour for the integral representation of delta functions. Therefore, one can also perform ensemble averages (| 1 . 26|) of 
trace relations with delta-functions of symmetry operators in their representation with coherent state path integrals, very 
similar to (11.4111. 14]) . However, it is important to distinguish between one-particle (11.1911.201) and two-particle operators 
(|1. 1711.271) in the various transformations to a nonlinear sigma model with inclusion of a coset decomposition for a SSB 
with a HST 
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2 Ensemble averages in model I and II for a normal-ordered Hamiltonian 

2.1 Precise time steps with shifts 'a£ p ' of the complex conjugated fields '^(tp+Atp) ' 

The ensemble averages of random potentials involve an additional, mathematical aspect which concerns the combination 
of fields of the two different branches (p — ±) of the time contour. This formal aspect is also implied by the original 
formulation of Anderson localization where one examines a 'return' probability of a test wave-packet (as an initial dclta- 
spikc at a particular space point). This 'return' probability measures the 'forward' propagation of the wave-packet, 



2.1 Precise time steps with shifts 'At p ' of the complex conjugated fields 'ip%(t p +Atp) ' 



away from an initial space point, and the corresponding 'backward' propagation of remaining wave-packet parts, back 
again to the same initial space point. Localization of wave-packets within a disordered potential is specified by a 'finite' 
'return' probability which approaches vanishing values towards derealization. According to this physical picture, we 
have to combine the doubling of bosonic coherent state fields ip$(t+) = ip3^(t), ipg{t-) = ijjg^(t), apart from their 
usual anomalous or 'Nambu' doubling (ipx(t p ) ■,4 ! *s{tp) ) of ordered systems with solely hermitian operators, with the two 
distinct branches of contour time. Consequently, there occur two different kinds of 'Nambu' doubling which we term 
'anomalous-doubled ordering' (|2.1ti2.4|) and 'contour time ordering' (l2.5H2.8p . the latter being marked by a bar under the 

doubled fields \£g(t p ) CUD, ^f iU) (EH), ^iSp) S^Tfe) (EH)- The 'anomalous-doubled' ordering (j2.H12.4p groups 

the fields ^§(t p ) = (iA?, P =±(t) (a = 1) ; tpg p=± (t) (a = 2) ) T according to complex conjugation regardless of the branches 
of contour time t p= ± whereas the 'contour time ordering' (j2.5H2.8p comprises fields of identical contour metric sign r\ p , 
regardless of complex conjugation. Aside from the 'equal time' and the hermitian conjugation of 'equal time' 'Nambu' 
doubled fields *§(t p ), *£°(*p) (j2.H2.2p or ^(t p ), *t' Q (t p ) (j2.5l2.6p . one has also to introduce time shifted versions of 

anomalous doubled fields ^>g(t p ), ^>h a (t p ) (|2.3l2.4p . &g(t p ), \t p ) (|2.7I2.8|) for the exact proper sequence of time steps 
within the considered 'quantum' problem which is represented by the coherent state path integrals of normal ordered 
Hamiltonians. The time shifted, doubled fields (|2.3l2.4p . (|2.7l2.8p are marked by the symbol ' ' to be distinguished from 
the equal-time, doubled fields (|2.1I2.2I) . (|2.5I2.6I) . Note that the underbar ' ' always hints to 'contour time ordering' 
(j2.5H2.8l) as in doubled fields, self-energy or source matrices and fields, independent of additional 'equal time' or 'time 
shifted' doubling 

'anomalous-doubled ordering' : t p = At 1 £ F > , Tq — a£ 1+ £ p ; 

(1) : ' equal time', anomalous-doubled field : 

^ (=1/2) (t P( =±)) = UsAt) , to,-® ; . **,-(t)) T ; (2- 1 ) 
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(1) : 'time shifted' At p , anomalous-doubled field denoted by ' ' : 
^ (=1/2) (i p(=±) ) = ( fa,+(t) , V^,-(t) ; fl+jt + At ) . ^-(t - At) ) T ; (2.3) 

a=l a=2 

(2) : 'hermitian-conjugation' with Hime shift correction' At p in the complex part : 

T,a 
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'contour time ordering' : t p = At 1 ~ T,P , , To — At 1+ £ p ; 

(1) : 'equal time', contour time doubled field : 

P =+ p=— 
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contour time ordering 

(2) : 'hermitian-conjugation' '^' of ^ equal time\ contour time doubled field : 

p=+ p=- 



£- Q( - 1/2) (t p(=±) ) = ( V4,+(t) , ^,+(t) ; #,_(t) , #,-(t) ) ; (2.6) 



contour time ordering 
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2 ENSEMBLE AVERAGES IN MODEL I AND II FOR A NORMAL-ORDERED HAMILTONIAN 



(1) 

¥ a( = l/2) / 
** (*p<=±)) 



(2) 

" a(=l/2) , 

£ 2 (V=±)) 
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contour time ordering 

'hermitian-conjugation' with Hime shift correction' 1 At p in the complex part 

p=+ p=- 
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(2.8) 



contour time ordering 

The Gaussian distributions (| 1 .411 .5[) with random potentials Vi(x), Vn(x,t) lead to the averaged coherent state path 
integrals Zj[3], Zu[3} which are composed of the same functional part (j2.11j) with the doubled source fields J^.g{t p ) 
(|2.12p and source matrices J%,.j{t P ) (|2. 1312. 141) (At p = r\ v At, At q = r\ q At, At > 0) 
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'ordering for anomalous-doubled source field' 
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2.1 Precise time steps with shifts 'At p ' of the complex conjugated fields 'ip%(t p +Atp) ' 



The one-particle part of (12.111) is listed in relations (|2.15H2.l7)) with Kronecker-deltas of time for precise time steps 
following from properly normal ordered Hamiltonians T and 'II'. We also give a 'laxed' kind of the one-particle part 
in relations (|2.16|2.17[) which is appropriate for a classical approximation, as a hrst order variation of fields, within the 
exponentials of the path integrals, but fails to result into the correct quantum expressions if integrations of coherent state 
fields and their complex conjugates have to be performed 
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(2.16) 
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According to the anomalous doubling, one can attain order parameter matrices (|2.18H2.2~lj) which are nonlocal in time for 
case T of static disorder and which are local in time for case 'IF of dynamic disorder. This corresponds to the described 
picture of a 'return' probability in Anderson localization for case T without an additional white-noise distribution in time. 
We can use the given order parameter matrices (j2.18H2.2T1) with their various forms for the anomalous doubled self-energy 
with 'equal time' doubling of fields (|2.18l2.20j) and for the 'Nambu'-terms extended density matrices with 'time shifted' 
doubling of fields (12.1912.211) . However, we can simplify to the spatially local case due to the contact disorder of spatial 
Kronecker-deltas with the second moments of the two Gaussian distributions (|1.4I1.5|) 
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anomalous doubled self-energy of static disorder ; (2-18) 

density matrices with anomalous terms for static disorder ; (2.19) 

anomalous doubled self-energy of dynamic disorder ; (2.20) 

density matrices with anomalous terms for dynamic disorder . (2.21) 



In order to achieve the precise, subsequent time steps at the integration boundaries, we have to introduce a modified time 
contour integration 'C", 'fg dt p . . .' which is extended to the parameter end points '0 P = —At 1+ ^ p ' and "T p — To+At 1 ~£ F '■ 
We thus take integrations with the extended doubled coherent state fields which are restricted to a single field (a = 1) at 
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which have to be combined with anomalous doubled matrices Mj?? g(t' q ,t p ) (|2.26l) for the proper exact sequence of time 
steps in order to obtain the correct propagation with the ensemble averaged generating functions (|2. 9H2.il]) . We list 
the various notations (|2.26I2.27[) of extended contour time integrations for the presented cases of disorder which are 
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straightforwardly generalized from the solely hermitian cases of interactions with the already described precise time steps 
and intervals in one of our earlier references 



.,6=1,2 



[dt p dt> q J2 E **S' q )Ml? t S' q ,tp)* a S P )= [ dtdt'J2 E % M£ rigiP tf,t) T) P %,(«); (2.26) 

C x,x' a,b=l,2 x,x' p,q=± 

/P+f a,b=l,2 
dt p dt' q j2 E *l'(*g)4w(C*p)*^)= m**'E E *^^ 9 (0^4;;^ P ( ^ ^ < )^np( ^ ) ( 2 - 27 ) 
x,x'a,fc=l,2 x,x' p,q=± 

= fjtpH E n(t P ) ■ 

• >C x a,b=l,2 

Concerning presentation and notation, the case 'IF of dynamic disorder turns out to have simpler relations than the case 
T because the static disorder involves the combination of two distinct time parameters t, t' in the self-energy matrices as 
in the already given order parameter matrices ^^"J *') (|2. 18[) and $jg ^ (t, i') (I2.19[) . The dynamic disorder with a 
white-noise time distribution restricts to a single time parameter in the self-energy matrices, as in the order parameter 
matrices &g p .g b q (t) (I2.20j) . ^g^.g'qif) (I2.21j) . but still incorporates the two metric sign labels p, q — ± as additional field 
degree of freedom for a 'disorder' quasi-particle. However, both cases T, 'IF of disorder allow for analogous treatment 
of HST's and coset decompositions of nonlinear sigma models, as one only considers the relevant reduction to stationary 
time states of fields with a single frequency 'ui p in case 'F of static disorder. 

Solely stationary states of time in case 'F result in the analogous transformations and derivations to nonlinear sigma 
models of case 'IF, if one replaces the single time parameter H' in ^ p .~ b q {t) (I2.20[) , ^ p .g b q (t) (|2.21l) with a single 
frequency parameter 'oj' in case 'F of static disorder. After the corresponding Fourier transformation to a frequency (or 

energy) contour according to (|2.28I2.29|) . we achieve for (|2.9[) the generating function Zj[3] (|2.30[) which contains two 
independent, non-stationary energy contour integrations 'fej p ' and ^huj' q in the ensemble averaged Gaussian part. As we 
reduce to stationary states or to a single frequency parameter lo = a/, we introduce the approximation (j2.31[) with various 

contour labels p, q very similar as in Zn [3] (|2.10[) 

Qmax = — ; < t p < +Tq ; To = T max ; (2.28) 

1 2ir 

< LU p < +ft Q ; il Q = Clmax = — ; Aw= — ; N t = T max /At ; 

AC J- max 



du p _ f +n ° dui+ f° du- 

ALU 



I 2-K \ 1 I 2.-K \ 

n ° diO + f+Vo+Au) duJ _ f +Cl P duJp 



(2.29) 



( 2ir \ ■ ■ ■ / ( 2-k \ ■ ■ ■ 2.^1 L i 2-k \ 

-Aw \tZZ> Ja ^T^> p=± J P \ T maa , ) 



0+ = -acj ; fl + = Q ; 0_ = ; £1_ = fi + Aw ; 
Zi[d] = (f^^J^J^] x (2.30) 

x exp f ?( hj^) rdUJp) e ^ p Atp Mujp) ) ( Scj^h rM) £l< Atq M<) 



Fbf>*,i>b 



Zi[3] ^ (FW^J^Jwitil x (2-31) 

The 'time-shift' correction of complex conjugated fields implies additional phases e 1 ^' A *p in the static disorder case whose 
consideration also leads to the exact, proper frequency or energy steps within the propagation of the time development 
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operators as in case 'II' of dynamic disorder. It is even possible to substitute the parameter R 2 H , according to relation 
(|2.32[) . and the contour time integrals by contour frequency integrals fg dui p / '(2tt /T max ) . . . (|2.29l) so that the nonlinear 
sigma model of case 'II' straightforwardly generalizes to case T of static disorder 

r>2 „. D 2 J *t 1L max /„ QrA 

K II ~ H I Jf ■ 

In the following sections [3J U] we can therefore concentrate on the dynamic disorder case with a single time variable in the 
self-energy matrix and can then transfer the result of case 'II' to the stationary case with a single frequency variable 'w' 
in the self-energy matrices for case T of static disorder. The approximation to a single frequency lo p for static disorder 
can also be attained at very later steps of transformations to a nonlinear sigma model, e.g. as one simplifies the two-time 
or two-frequency dependent, nonlocal disorder self-energy matrix to a single dependence with contour frequency tu p . 



3 HST for dynamic disorder and repulsive interaction with 'hinge' fields 
3.1 Anomalous doubling of the one-particle part 

We perform the anomalous doubling (|3.1|) of the bosonic fields with inclusion of the contour time metric rj p , due to the 
ensemble average with a dynamic disorder. This defines a 'Nambu' metric tensor K pq (|3. 21) for the anomalous doubled 

fields ^g b q {t), ^3 p {t) whose dyadic product determines the density matrices R^ b pq {t) (|3.3|) with anomalous extended parts 
in the off-diagonal blocks a ^ b. The quartic disorder interaction of bosonic fields is therefore equivalent to the trace 
relation (13.41) with the density matrices R~ b pq (t) Q3.3P which are modified by the 'Nambu' metric tensor p.2[) with a final 
multiplication of a factor ' j' corresponding to two dyadic product operations of (|3.1|) with anomalous doubling of relation 



rSp+^p) v P Mt P ) = \ (rsp+^p) v P Mtp) + Mt P ) vp rSv+*tp)) = \ (*?) Kl nth) ; (3.1) 

b 

= 5 a b S pq T) p = S (l b S pq diagj +1 , -1 ; +1 , -1 j ; (3.2) 



K ab = ( 

pa I / ~\ 22 



'PI / PI 0=1 a=2 



(^PJ tq) 

( 



\ 
I 



^ lfe,-(t) J 

r s +(t+At) 

(*- a *) 



1a=l 

v 

a— 2 



/ 



6=1 



V>|,+(*+Ai) ^,+(*+Ai) VI,+(^+At) ^,_(^-At) 

^,_(t-Ai) ^i + (t+At) ^_(«-Ai) ^,-(i-Ai) 



21 



fc=2 



(3.3) 



Jt-At) i> S ,+(t) _(*~At) ^,-(*) 



22 



/«// 



P9 ■ 



- Tr 

4 a,b;p,g 



R% q {t) K R% qp {t) K 



(3.4) 



This anomalous doubling has also to be taken for the one-particle part (|2.15[) so that we accomplish relation (|3.5[) with 



the anomalous doubled fields _(i')> ^£ p (i) and the doubled one-particle operator 5C|? q -^ v {t',t) P-61) whose lower 
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block diagonal part - (if, t) (I3.8[) follows by transposition from the upper block diagonal part "Kgi - (t', t) (|3.7|) . 

Corresponding to propagation with the exact, precise time steps, we outline these two one-particle parts in Eqs. (I3.7l3.8f) 
with the exact Kronecker deltas of time and also describe a 'laxed' kind of one-particle operators in (|3.9|3.10|) . only 
applicable in classical approximations or equations 

[ dtpj^^iitp + tip) Hp{x,t p ) i/j s (t p ) = J dtpY^Mtp) Hp&h) *P*s(tp + tip)= (3.5) 

^ X X 

= Jjt p dt' q J2~¥*(t' q ) K b s % 3 (if q ,t p ) n(t p ) ; 

^ x,x' 

&l' ; a(*£> t p ) = diag (jfc^. 2 (i£, tp) ; &&.g(t' q , t p fj ; (3.6) 
X$. s (t q ,t p ) = ( - ih g(^zfaz^jz£giz^ + (h(x') - t g p ) S(t' q - ^ V t?, fo >g ; (3.7) 

&I^,*p) = (^;x(*' 3 '*p)) T = ( -ih 6 ^-"'^-^-^ + {h T (x") - i e p ) 8(t' q - tp^S qp Vq Svj ; (3.8) 



^, s (t' q ,t p ) ~ W|^(t' 3 ,i p ) = diag H+(x',t' + ) , ; ,*' + ) , -HH#,1L) )S SP)S S qp S(t' q - t p )(3.9) 



i? p (x , tp) = l£p+ 2m + U ^ ~ ^° ' H p(%, tp) = + lh Ql * £p + 2m + U ^ ~ Mo ' ^ 3 ' 10 ^ 

3.2 Anomalous doubled self-energies and their coset decomposition 

In the following we convey the results of our earlier articles to the presented case 'II' of dynamic disorder with the 

fab 
P9' 



additional contour time metric tensor rj p and have to use a modified 'Nambu' metric tensor K pq , K ab (|3.11I3.12[) which 
changes the anomalous parts in the self-energy to anti-hermitian relations 

K ab = [ IU W 22 j ^ Sab 6pq 7]p K ab = S ab S pq diag{+l , -1 ; -1 , +l} ; (3.11) 



^ 'PI 7 PI a=l a=2 



n ab = S ab diag{^; ^} (3.12) 



a—1 a— 2 

M 



The final HST is achieved by one half with the self-energy density variable cr R ' (x, t) p.!3p and by one half with the 

'Nambu' parts St^ pq (t), St%} pq (t) (a ^ b) P~T3) of the anomalous doubled self-energy matrix &tf pq (t) (|3~T1) whose 

block diagonal, hermitian self-energy density parts <5E?j! pq (t), 5Y?~ pq (t) p.!5|) are related by transposition and are only 

used as 'hinge' fields in a SSB with a coset decomposition. The off-diagonal, 'Nambu' matrix parts SEg 2 pq (t) 7 <5£|? _ (i) 

p. 161) of S^i b pq (t) p.!4p are symmetric matrices and are related by hermitian conjugation. One thus acquires following 

real and complex parameters SBg- :++ (t), 5B^ (t) and 5B^ (t), SBg. |_(t) = 8Bt . , (t) and the solely complex field 

variables Scs-++(t), Scg. (t), Scg- |_(t) = Scg.^ (t) in the corresponding anomalous doubled part (|3.16p of the total 

self-energy (|3.14[) 

a%l(x,t) e R ; 5t\) pq {t) ■ 5±% pq (t) ; dt% q (t) ; S±l) pq (t) in analogy to R% q (t) ; (3.13) 
5± ab (t) = ( S ^pi^ 5 ¥%pi® | • (3 14) 

x -' pq{ ) { sti^® s^ pq (t) ) ' {3AA) 

syii m f SB s , ++ (t) SB s , + ^(t)\ ~ 22 / SB Sl++ (t) 5Bi + _(t)\, 

^x M W - ^ SBi + _(t) 8B 3 .— (t) ) ' fl2j *3*W - ^ 5B s .+_(t) 8B s ,__{t) ) ' {6Ab) 

St%£(t) = Stl] pq (t) ; SB s , ++ (t) , SB S ,__ eR ; SB S]+ _(t) eC ; SB £ ,_ + (t) = SBl + _(t) ; 



3.2 Anomalous doubled self-energies and their coset decomposition 
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Ay 12 

02 ^S;pq 



(t) 



it) 



-(*) 
-(*) 



5Ef* (i) = 5£^+(i) 



ryl2,T _ ,A12 



(*) 



AY 21 

0Zj x;pq 



e C 



(*) 



*£;+-(*) 



^4 

fcj? ;+ _(i) = 5c £; _ + (t) . 



.(*) 



(3.16) 



We combine the self-energy density variable o^ u {S,£) (|3. 131) and the anomalous doubled self-energy matrix ST,~ b pq (t) 
([3711} into (|3~T71) with block diagonal parts E|f 9 (i) f[37T8|) . £§? pg (t) (|3~T9l> and introduce the modified self-energy matrix 
S^g b pq (t) K (|3.20I3.21|) with 'Nambu' metric tensor 'K' (|3.11l3.12p and anti-hermitian related, off-diagonal, anomalous 
doubled blocks St^ q {t) = i 6£ a s f pq (t). This allows to perfo rm a coset decomposition into densities and bosonic parts 
according to Sp(4) Sp(4)/U(2) '<g> U(2) 



Si 1 

x;pq 

SY 21 

x\pq 



(*) 




(*) = 



= a£(x,t) X + 

\ w ^x\pq 

= -VRu&t) % 5 P<I + <^t;pg(*) 

/ f SB s , ++ (t) -SB s . + _(t) 

^ ; ++(*) + 



<5£y 



(*) 



SY a -f b (t) = i 6± a -? b (t) ; 

/ x;pq\ > ' 
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-5Cjf ; -H 



(*) 



8c 3 . + -(t) 
Scg- — (t) 

5Bi + _(t) 
SB s .^(t) 
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/<</ 
2:2 



(3.17) 

(3.18) 
(3.19) 
(3.20) 

(3.21) 



As one includes a further, derived 'Nambu' metric tensor I pq (|3. 221) so that the various kinds of self-energy terms with 
hermitian or anti-hermitian off-diagonal block parts can be transformed into each other p.23|) . one finally succeeds into 
the coset decomposition (|3.24[) with block diagonal densities <5£^ pg (x, t) K (similar to (|3.15|) ) and coset matrices T pq (x, t) 

(|3.25[) . The latter coset matrices (|3.25p consist of the generator Y pq (x, t) (|3.26[) with sub-generators X pq (x, t), —X pq (x, t) 
p. 271) as the 'Nambu' or anomalous doubled field degrees of freedom 



fab 
1 pq 



Sab S pq diagj^ 



-1, +1 



-l , 



:} ■ 



(3.22) 



f ( ^MpqW S ^s 2 , P q^ I f = a ( - 0) (xt) K + SY a - b (t) 



f (f , t) = exp { - f p f (f , t) } ; f - 1 (f , t) = cxp { f p f (x, t) } ; 

V)" (X m {x,t)f 2 \ ( (o) 



y_ (a;, t) 



X pq {x, i) — 



li 



pi 



X*(x,t) 



i] p X\ (x,t) n q 



(x pq (x,t) 

\ / na 



12 



-5czj ;++ (a?,i) (5cd ;+ _ (x , t) 
-<5co ;+ _ (f, t) 8c D . — (x,t) 



; -i] p Xl q (x,t) n q = 



6c *D-+- 



Sc 



D-A 



_(x,t) -5c* D . + _(x,t) 
_{x,t) -5c* D . (x,t) 



(3.23) 
(3.24) 

(3.25) 
(3.26) 

(3.27) 



It is further possible to diagonalize the various block diagonal density parts 5T/]j. pq {x, t) K as in (|3.28[) to (|3.35p and 
the various 'Nambu' generators Y pq (x, t), X pq (x,t) as in (13.36)) to (|3.40p . In the case of the density parts we have the 
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diagonal eigenvalue elements (|3.29[) with diagonalizing 'rotation' matrices (|3.30p so that the lower block diagonal '22' part 
is related by transposition to the upper '11' part 



5V£. m {x,t) K = Q; p Y a (x,t) SA a p ,(x,t) Q a p ? q (x,t) 

( ( SB D]++ (x,t) -SB D]+ -(x,t) x M 
SB* D . + _(x,t) -5B D .—(x,t) 



SA a p (x,t) 



Qpq 



^S-pgO^J t) K qq 



V 



p-8X p [x,t) 



\ 



-5B D , ++ (x,t) 8B* D . + _(x,t) 
-5B D .+-(x,t) 5B D; ^(x,t) ) m ) 

—p-5\p(x,t) 



diagj (+SX + (x,t) , -6X-(x,t)) ; (~SX + (x,t) , +<5A_ (£,*)) \ = 5X p (x,t) Kj 



^A p (x,t) = diag{£A + (x, t) , 5X~(x, t)} ; SX p (x, t) £ M ; 

= (exp{«£ D (z,t)}) ; Qf q (x,t) = (exp{iS5(i,t)}) ; 







-S2, ;+ _(j?,t) 



exp{i/3£)(f,t) } ; ®_d ;+ _(ie, t) <E 



Ul.-l/ 



(x , i) 



Q^(2, t) 



O 22 '- 1 ^,*) (-p-SAp^.t)) Q 22 (x,i); 

J D;pq\ 



6t^(x,t). 



(3.28) 



(3.29) 

(3.30) 
(3.31) 

(3.32) 

(3.33) 
(3.34) 
(3.35) 



A similar diagonalization (|3.36I3.37|) is achieved for the off-diagonal block parts with sub-eigenvalue parts Xo- tPq {x, t) (|3 . 38[) 
and 'eigenvector' matrices P^^x^t), P^^i^^) P-39j) in such a manner that the symmetry relations of (|3.25p to (|3.27l) 

are still retained. This involves the complex parameters c++(x, t), c (x, t) as eigenvalues and Cdh (x,t), CJ>j (x, t) 

as the angular parameters of the rotation for the off-diagonal elements 5cd,-\ [x, t), 6c* D .^ (x,t) and diagonal elements 

5c D . ++ (x,t), Sc* D . ++ (x,t), Sc D . — (x,t), Sc* D . (x,t) within X pq (x,t), —r) p X^x^) r\ q (I3.27|) 



Y pq (x, t) 
Y D . ixi (x,t) 

X D - pq {x,t) 
c ++ (x,t) 
^2x2 

GD;2x2(x,t) 



P ixi (x,t) Y D . Axi (x,t) P ixi (x,t) ; 
11 

I 

' pq 



( (oY X Dm (x,t) 

\ / pq 

{ -Xl pq (x,t) (V) 



22 

PQ 



Pix4,(x, t) 



Pi x2 (x,t) 







^2 x 2(^1 1) 



-c ++ (x,t) 
c—(x,t) 



~xl>. tPq (x,t) — ( 



c* ++ (x,t) 

-c*__(x,t) 



pq \ \ ' / pq 

\c+(x,t)\ exp{i (f+(x, t)} ; c (x, t) := |c_(x, i)| exp{z<^_(£", t)} ; 

22 



CXp|« 6^:2x2 (2,t)| ; P 2 2 X2(^*) = GX P{* C5;2X2(^0} ; 



GD ]Pq {x,t) - ^ ej,. + _(f,t) 
16^(2,*)! exp{i7 D (x,i)} . 



-e 0;+ _(x,t) 



Ax 2 C 2 -)*) — P2x2( x ,t) 



l"l 



(3.36) 
(3.37) 

(3.38) 

(3.39) 
(3.40) 



3.3 HST transformations with 'hinge' -fields 
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3.3 HST transformations with 'hinge '-fields 

Eventually, we can collect the various parameters 5Bg- ++ (t), SBg ; --(t), SBg-^ (t), SBi^ (t) for the density parts 

(|3.43|3.44|) and Sc 3 - pq (t), 5c* s . pq (t) within the 'Nambu' terms ([3^45]) in order to determine the HST (f3T46|) 



4lm 



'hinge' functions : SE^Jt), 6T,f (t) 



5B S , 



SBsk. 



SBl + _(t) 



ry22 



(*) 



+■ 

+ (*) 
(*) 

+ (*) 



SB S 
SB s . + _(t) 

5cg ;++ (t) 

6c S ; + -(t) 



SBg, + _(t) 
6B s .__(t) 

8Bt + _(t) 
SB s .__(t) 

6c£ i+ -(t) 
Sc s - — (t) 



6B s . + _(t) e 



-.22 



(0 



(3.41) 
(3.42) 

(3.43) 
(3.44) 
(3.45) 



The HST of a quartic interaction of fields can be taken in various manners; we consider the case where one half of the 
quartic interaction (with correspondingly reduced pre-factor R]j/(4h 2 ) in the exponent) is transformed by the real self- 
energy density variable cr^j (x, t) as invariant vacuum or ground state in a SSB and where the other half of the quartic 
interaction (also with pre-factor i?| J /(4fi 2 ) in the exponent) remains within the anomalous or off-diagonal blocks of the 
Gaussian transformations. This can be verified as one decomposes the trace inside the exponential of the last line of 
(|3.46[) . As we remove the trace and dyadic product inside the exponent of the last line in (|3.46j) . one notes that terms 
with the block diagonal self-energy density 5Z3^J (t), ^ s 2 pq (t) cancel and only the self-energy variable cr^^x^t) and 
anomalous parts Sf^ s 2 pq (t), ST,g} pq (t) couple to the anomalous doubled, bilinear fields $?g b q (t) . . . ^%. p (t) 



(x,t)] exp | - \ ±- £ T dt £ 4°i & *) 4°i % t) 



(3.46) 



d\E$L. (t)K] exp 



x exp ■ 



i 



- / dt 



1 1 

8 Rjj 



Tr 



dt 



V Tr 

_ a,b;p,q 



SZ±Jt) K K 




-nut) 




K 



In analogy the HST is obtained for the repulsive, quartic contact interaction of Bose fields; one half of the HST follows 
from a Gaussian identity with the real self-energy variable ay\x, t p ) and the other half is given by a Gaussian identity of 

self-energy matrices Sai b (t p ), 5<jg h (t p ) which are shifted by the anomalous doubled density matrix Rg b pp (t). We emphasize 
the missing couplings of the different branches '±' of the time contour, due to the hermitian property of the repulsive 
contact interaction, so that one has to adapt the corresponding symmetries of Sag (tp), Sai b (t p ) to that of a single time 
branch with 5ai b (t + ) and Sai b (t^) being unrelated self-energy field variables 



exp 



d^voix^p)] exp 



x / d[5af(t p ) k] exp 



b 

(t p + At p )^j (lpg(t p ) 

L dt 'V°- 

X 

— I dt p y~* tr 

4ft Jc ^ a,b 



(3.47) 



i 

2ft 



'V ( x >tp) cr V ( :E '^p) 

Vo-isp 
Saf(t p )Kaf(t p )n 



Vq 
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°"I 2 (*p) 

8af{t p ) 



x 




tr 



R^ppi*) 



5af{t p ), Saf(t p 




°¥(t P ) 5ag(t p ) 
Svfitp) -°f{t P ) 



0y o OMp) G 



S5g(t p )=t6vg(t p ) (a^b); 
8ol\t p )=5of{t p ) . 



(3.48) 
(3.49) 
(3.50) 
(3.51) 



3.4 Removal of the 'hinge'-fields from the generating function 



As we insert the two HST's (j3.46|3.47p of the quartic, non-hermitian disorder term and the repulsive contact interaction 



into (|2. 1012. 1112. 311) , we achieve the ensemble averaged path integral Zn [8] (j3.52[) with only linear and bilinear anomalous 
doubled fields q (tf), &g p (t). Aside from the linear coupling to the doubled source fields J^. s (t p ), Jl^(t P ), we 
abbreviate the bilinear term of fields ^t b , it') ... tyt (i) by introducing the matrix M|? s (t' q , t p ) (|3.53p which consists of 
the one-particle part Hgi^(t' q ,t p ) (I3.5ll3.10p . the source matrix g(t q ,t p ) for generating bilinear observables of bosonic 
fields, the condensate seed field J^p.^(t p ) (|2.13l2.14l2.27p and the various self-energy variables cr^(x, t), <jy^(x,t p ) and 
matrices <5£| a (t), Sai a (t p ), defined in previous sections 



Z H [3] = / d[S^ pq (t) K] exp 



1 1 

8 R 2 H j„ 



tr 



+T 



dtJ2 Tr \s^ pq (t) K 6^ qp (t) K 



8at{t p )n6df{t p )K 



Vq - ie p 

Vq—IS v 



x J d[S5f {tp) k] exp | ± jf dtp £ t 

x J d[(T^(S,t p )] expiry dtp ^ 

x | d^, (f , t)] exp | - \ ±- £ T dt £ a£ (f, *) a W (£ , t) | 

x / d[^(tp)] exp{ - ± / +T ^ df J] g (0 ^.^(f.t) 



x exp 



i 

2fi 



H ^pE( J fe(*p) *t(*p) + **(*p) ^(*p))|; 



8 b £>s{tqitp) 



M b ^ s (f q ,t p ) = 0<ll s (f q ,t p )+K x ' x ^ if+ 

+ *(i - O S pq i lp hl%, s (t p ) + S ab ^ tp) + i ff w (fj ^ 



ffi^(t) + 2 5p 9 77p fci 1 (t p ) Sm p (t) + 2 5p 9 Vp 5af{t p ) 
+ 2 6 pq Vp 5af{t p ) - (5±l% p (t) + 2 6„ Vp 5af{t p ) 



(3.52) 



(3.53) 



K 



After integration over the bilinear, anomalous doubled fields in (|3.52l) . one attains the inverse square root of the determi- 
nant of M|» g{p' q , t p ), which is transformed to a '— ^tracelnf. . .]' in an exponential with normalizing term H£l max J\f x and 



3.4 Removal of the 'hinge' -fields from the generating function 
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to a propagator part M-, 1 ^ (t', t) between bilinear, 'Nambu' doubled source fields J]^.^,{t' q ) ■ ■ ■ J$-${tp) 



Z n [3} = I d[6Z£ tPq (t) K] exp \ ~ \-L I dt^ Tr 6V$ pq (t) K 8Z% qp (t) K 



8R 



II Jo 



J d[5df{t p ) k] exp | A J dt p £ tr 
d[o^ (x, t p )] exp | ^ J dt p 



_ a,b;p 7 q 
x 

Sai b (t p )HSa b s a (t p )K 



Vo—lEp 

V -ie p 



(3.54) 



| rf^ (x , t)] exp | - \ ±- jf dt£ 4°i (f , t) CT £ (*, t) 
exp { - i £ ^r/ p £ m max jV x Tr In [m|? 5 (t' 9 , t p ) 



a,b\p,q 



x exp 



2 h 



" s*' p,<?=± 



In order to simplify the self-energy matrices in a coset decomposition, we shift the matrix SYfi pq (t) and self-energy variable 

Oy^(x, t p ), according to (|3.55|3.56j) and p.57|) . respectively, and perform the transformations p.58|) of t' q ) (j3.53j) 

with the 'Nambu' metric tensors p. 213.1113. 1213. 221) to a modified matrix Mf g ,(t p ,t q ) (I3.59l3.60j) . This does not alter 
the value of the determinant (|3.61[) and allows for the factorization (j3.63H3.66j) into density-related self-energy matrices 
and coset matrices T£ q (x,t) (j3.25H3.27j) . Similar transformations (j3.62j) are also accomplished for the propagator part 



from M^r(t p ,f q ) to MZ^ D (t p ,t' q ) with 5Y,f pq (t) K (13331 for the coset decomposition Sp(4)/U(2) ® U(2) 

5tt tPq (t) -+ 8t% q (t) - 2 5 pq r) p 5Sf (t p ) ; 

^S;pq(t) ~> ^g;pq(fy — 2 S pq T] p J^.g(t p ) ] 



^x,x' (tpi t q ) 



&x,x'(tpjtq) ~ •^■x*,x' ftpi tq) K 

DET{M^,(i pi ^)} 



k r 1 [I K M a s %(t p ,t' q ) K i) r 1 K 

k r 1 (i k Mf t g,{tp,t' q ) k i k\ k r 1 k ■ 
v — — „ ' 

i k Mf s ,(t p ,t' q ) k i k 

Hp(t p ) + a$(x, t p )\ S ab S pq 5{t p - t' q ) S s , s , + 

6±% pq (t) i6±% pq (t 
^6i% pq (t) 5±% q {t) 

= H p (t p ) 5 ab S pq S(t p - t' q ) 5 3 ^i ; H p ~ 2 (t p ) — ^i?p _1 (i p )^ ; 

= DET{i: r 1 Mf iS ,(t p , t' q ) K t 1 k) = DEt{m^,(* p , t' q )) ■ 

= k ikM^j? b (tp,t' q )ik; 



ik+\8{t-t')8^. 



^pqW *St%„{t 
l5±f(t) 8tf(t) 



ah 



K = 



f$(x,t) Stir/^t) K t£f\2,t) 



(3.55) 
(3.56) 

(3.57) 

(3.58) 



(3.59) 
(3.60) 



(3.61) 
(3.62) 

(3.63) 
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3 HST FOR DYNAMIC DISORDER AND REPULSIVE INTERACTION WITH 'HINGE' FIELDS 



-8%D;pq&t) T] q 



5BD-++(x,t) -5BD-,+-(x,t) 
5B* D . + _(x,t) -8B D -—{x,t) 

-SB D . ++ {x,t) 5B* D . + _{x,t) 
-6B D]+ _(x,t) SB D '.,—(x,t) 

(0), 



(3.64) 

(3.65) 
(3.66) 



Since the shift of the matrix 8T 1 '§'. pq (t) (13.551) and self-energy variable <jyj (x, t p ) (|3.57|) has removed the self-energy matrix 

Saf(t p ) and variable a ( ^ u {x,t) from Mf 3 ,(t p ,t' q ) or Mf g ,(t p ,t q ) (|3.59I3.60|) (compare with f3~5^) ). we can completely 

eliminate latter self-energies by integration of Gaussian identities which involves the appearance of a new parameter /Jp 
(|3 . 69[) with the ratio Rjj/(hVo) in the denominator 



dWnl (3 *)] Cxp { ~ J 6 d< E a Rn ( f ' *) a[ Ru ^ ' *) j Gxp { ifi J G dt P E 



+T 



exp 



1 



II Jo 

1 p+f 



Vq - I £ p 



Vo-lE p 



exp 



1 1 

8 R]j 



V -i£ p 

1" dtT Tr (S^ pq (t) - 2 6 pq Vp Udt(t p ) + Jf^Sp)) ) K 

_ a,b,p,q |_ ^ V ' ' 



(3.67) 
(3.68) 



x [d^ qp (t) - 2 S pq Vp (§5 b / (t p ) + J^(i p )) ) X 



exp 



x exp 



x exp 



1 1 

8 R 2 U jo 

+f 



£ p,q— ± 



a.b 



2 -R// Jn 



2 -Rf j jo 



+ T 



W (7/) = 



i 



a . 6 



tr 

a.b 



;pp\ b ) n u inp 
J4np;x(tp) K Jip^gitp) K 



(l - (*/2) ?7 P (i?? / /(^o) 



(3.69) 



After insertion of the Gaussian identities (|3.67I3.68[) into the path integral fl3.54[) . one obtains Zjj[3] p.70p which is 
determined by the self-energy variable Oy^{x, t p ) and the self-energy matrix term ST,i b pq (t) K for the coset decomposition 
as remaining field degrees of freedom 



Z n [S\=exp 



2 R 2 U jo 



*E E f 1 - 4 m ) K Widtp) « 4%dt. 



x j d[a$(x,t p )} 



x p=± 
y2 



(3.70) 



3.4 Removal of the 'hinge' -fields from the generating function 
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x exp 



i 

2ti 



Jc 



V -lEp 



x / d[5V$ pq {t) K] x 



exp {"^ /^E ^ (l-^M^) tr [S^ pq (t) K 6^ qp (t) K 



x exp 



1 f dt 



^ Vp V nn max A4 tr In t p ) 

' ^ — ' a,b 



x exp < - 



ft 2 

max 



2 /i 



r dtdf ee^* e ^(4)^ M ^ ba (4^p) ij ^)} • 

p,«3=±x,x' a,b=l,2 J 



Finally, we factorize the self-energy matrix term S^g b pq (t) if inside M^ & -,(t p , t' q ) (|3.71|1 corresponding to the coset decom- 

and apply this factorization inside the action Adet[T, <5S_d, ct£ ; 3] (|3.73[) from 



position Sp(4)/U(2) O U(2) 

the determinant and within yi 7 t 7 ^ [T, cr^ 1 ; 3] (I3.74[) from the propagator part of the source fields J^.^(t p ), J^".g'(t q ) 



Mf^ = [i? p a (tp) + <xg(z,tp)J 6 ab S pq S(t p - t' q ) 5 St3 , + I us '*'^'"<i> ik+ (3.71) 



+ 



i 5(t - t') 6g t g> f(x, t) Q 4 ^ A {x, t) SA ix4 (x , t) (f(x, 1) Q ixA {x, t)^j _ 



T(x,t) Q 4 ^ 4 (f,t) = T (f,t) ; T Q x (x,t) = Qi X i{x,t) T 1 (x,t) 
r (°). 

■)2 /-+T 



5 £ H*^ E m rna X tr In [m 6 ,? ip) 



o,b=l,2 



(3.72) 
(3.73) 

(3.74) 



In subsequent steps we straightforwardly outline how to factorize the matrix M^ b -,(t p) ^) (|3.75I3.76|) with the coset 

matrices T pq (x,t) (|3. 22113. 27j) and with the various 'Nambu' metric tensors (|3. 213. 1113. 1213.22]) into the modified matrices 

^tA^A^^D K) (EZZD, NZK,{t p ,t' q -J- 1 5± D t 1 ) GEE), N£%,(t P ,t' q ) (023, 6 a s b s ,(t p ,t' q ) gSJJ where the first two 

kinds of matrices (|3.77l3.78p do still contain the block diagonal self-energy density matrices <5S^. (x, t), <5£ff. (x, t) as 
'hinge' fields in a SSB with a coset decomposition Sp(4)/U(2) ® U(2) 



= [f{x,t)\ t X a s %(t p ,,t' q ,;S± D K) (f-^x 1 ,0J , 



b'6 



(3.75) 



MgVM) - * * ^I'(V>^;^ if) KI-^K) KIK f^ b "(x',t f y, (3.76) 



m'*:(t P >,t', -J- 1 s£ D i- 1 ) 



Ni%,(t p , t' q , r 1 s± D i- 1 ) = (k r 1 N^,(t p , t' ql s± D k) k r 1 k 

NlAtp^I- 1 SZ D I' 1 ) = Nl b s ,(t p ,t' q ) + 5(t-t')5 s , 3 , Kl d ™ { ' ] 

m% (t P , t' q ) = k\( \k;(t p ) + a$ (x, tp )] s ab s pq s s , s , s(t p -t' q ) + 



ah 



-6±%<x',t') 



K 



(3.77) 
(3.78) 
(3.79) 
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3 HST FOR DYNAMIC DISORDER AND REPULSIVE INTERACTION WITH 'HINGE' FIELDS 



(f(d 



-l:ab' 



pq' 



x,t)l) ; H b q ,(t' q ,) (f(x',t')l) _ -H;{t p )8 ab 5 m 5^5(t p -t' q ) ) \{t p ,t' q ) + 



b'b 



q'q 



ab 



K < (T I 



ab 



1 iJLiK[fi)\(t P) t' g ) = K r l otAhA) J ; 

x.x' 



tp) = + < } <)J «6- V 5 (4 - *p) + (r- 1 Ijj-JKT^ {t' q ,t p ) + (3 

f£ i6 °V> ^'(V) r p % a (x, t) - <5 fca V 8{t' q - t p ) 



The factorization with the coset matrices T£ q (x, t) 



S'K^" 1 ,T) 



has isolated the 'hinge' fields for the density related terms 



and has introduced a new gradient term <Wf(T , T) for the anomalous doubled parts 

^{T' 1 ,T;t' q ,t p ) = [T~ p )' M (a?,t') (V) T p' P {^,t) ~ H q(t' q ) ha S qp 8 S ,. 3 5{t' q - t p ) 



(3.81) 



so that we can use this separation into 'hinge' density-related self-energy terms and 'Nambu' gradient terms for the actions 

r$;3] (£72} and A^ tJ jf,S±r 
(|3.79l3.80p without 'hinge' fields are related by 



•Adet[T, ££d, (Jy ; 3] (I3.73[) and Ajt [T, ST,d, Vy , 3] P-741) . Since the determinant and propagator of the matrices 



DET 



(N$; s ,{t p ,t' q ) 



DET 



N^ ab {t p ,t' q ) = i- l bzy>(t p ,f q )iK, 



we can considerably simplify the actions (|3.73I3.74I) by following transformations 



DET [M ab (t p ,t')) = DET(NZ%(t pi t'-r l 5Y, D r l ) 



A D ET[T,S'E D ,a^;3} = \ /, y % )^ ^mw -A4 



(3.82) 
(3.83) 

(3.84) 
(3.85) 



x tr In 

a, b 



N§ b s ,(t p ,t q ) + 6(t~t>)6 s , s ,K D --™ 



K 



M-tv ■ *i) = ( f ' *) X K ^ (V - £ ; ^ SZd I- 1 ) K I- 1 f^ b \x", f ) ; (3 

r+T o,5=l,2 , 

•'0 a p,9=± 99 



o 2 



2h 



(3.87) 



x * N£$ a (t^X'-J- 1 SZd I' 1 ) K (I' 1 t~\3,t) I) J%.S P ) ; 

r 1 f-\x,t) i) JIAtprf = & P ) O' 1 T-\x,t) /V = j+utv) (i k f(x,t) k r 1 



(J- 1 f- x {x,t) fj = r 1 expjy^x,*)} 



I = r 1 exp 



+fj X fi 

-xt o 



I = 



(3.88) 
(3.89) 



r 1 K exp 



X 
-fjX^fj 



K 1 = 1 K T(x,t) K r 1 



After substitution of above separating relations into self-energy densities <5S^ a p(J (x, t) and gradient term (|3.81|) of coset ma- 
trices into the path integral (|3.70|) . we can 're-introduce' integrals of bosonic, anomalous doubled fields ^^,(t' q ) . . . ^%{t p ) 
so that the determinant and propagator part are changed back to the bilinear, 'Nambu' doubled term of bosonic fields 



3.4 Removal of the 'hinge' -fields from the generating function 
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and to the linear coupling with the condensate seed fields J^,.^(t p ), J^ a $(tp)- br this manner one performs a projection 
onto the anomalous doubled field degrees of freedom 



Z n [d] =exp 



2 R 2 U jo 



+T 



x p— ± 



(3.90) 



x / d[o-^(x,ip)] exp 

x exp {^//^E 



1 Rj T 



d[<yEs'(t) ^] 



x exp 



x exp 



1 1 

8Rfi 



Vo-ts p 

a: p,g=± 



+T 



+ T 



x exp 



2 R 2 n jo 



I 
2ft 



^EE^f 1 



x p=± 



*E§L(t) « ^ iaf (tp) k x / d[Va(t P )] x 



tr 

a . 6 



, ay, 11 n 



KN^Jt'a, tp)K + 



^11 




Ay 22 

OZj D;2x2 



n bo 



(Op) *S(*P» >< 



x'x 



dt p dt' q j2^P(t' q ) 

x*,x f 

cxp {~4 /.^Ef^Scv) (iKf{x,t)ki- i ) a mtp)+*t a (t P ) i^- 1 T-\x,t) i) j* ;af (v) 



•^. >a (**) 



However, the presented factorization of 'M'gg, (ipjig) with the various 'Nambu' metric tensors and coset matrices has 
modified the original path integral (|3.70|3.71|) or similarly (|3.73|3.74j) to (|3.85I3.87|) or (|3.90|) in such a manner that the 
bilinear coupling of doubled, Bose fields ^%{t' q ) ■ ■ ■ &g(tp) with the self-energy densities SY^.^x, t), 5Y,'j^. pq (x, t) cancels 
in (13.90)) and restricts the corresponding exponential (|3.91[) to the value of unity 



exp { ~ in Jc dtp dt ' q E^ ^ 



Ay 11 

02 "D;2x2 



-6E 



22 

D-2x2 / x',x 



ha 



(3.91) 



The above transformations with 'Nambu' metric tensors and factorization of the total self-energy ST,t b m (t)K have resulted 
into a projection onto the anomalous doubled field degrees of freedom with the coset matrices T^{x, t) and the self-energy 
variable <jy (x, t p ) as the invariant vacuum or ground state in a SSB. We thus obtain the path integral (I3.92p with the 
operator 6|° g(t' q ,t p ) (|3.93[) . being composed of gradient 6$C(T~ 1 ,T) (|3.81|) and density part H p (t p ) + Uy^ (x,t p ), where 
the block diagonal 'hinge' densities SYljf. pq (x, t) are only contained within Gaussian factors and the traces of 5T,i b pq (t) K 
or ST,i b pp (t) k so that these 'hinge' fields can be removed by Gaussian like integrations 



Zii[3] = expi - 



2 R 2 H «/■ ) 



+ T 



— — v ' ' a,o 



x / d[a^(x,t p )] exp 

'{a/a*? 



p=± X 

1 R? n 



(3.92) 



x exp • 



x exp 



4(hV Q 



V - 1 s p 



d[SZf pq (t) K] 



1 1 

8 Rjj Jo 



dt E E I 1 - ^ $ 

p,q=± x 



II 



tr 
a . b 



8T,f(t) K ST%(t) K 
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x exp 



x exp 



x exp 



2R 2 n JO 



+T 



p=± X 



tr 



1 

~ 2 



2 ft 



^77 P V m max jV x tr In i p ) 

a; 

_l_f- a' 5 6'=l,2 a 5 6=l,2 

E E E^4l^)^ f ^^ 

p',q'=± p,q=± x,x' 



O b xlx(t' q ,t P ) = (H b q (t' q ) + a$(x , ,t' q )) 5 ba 6 qp 8&,g 6(t' q - t p ) + (T" 1 I I K f)^ Jt' q ,t p 
+ fc 1;bQ '(x', O #£(V) r p % Q (f, t) - H b q (t' q ) 5 ba S qp 6 3%s 6(t' q - tp) 



3 



(3.93) 



5M(T- l ,r) 



3.5 The change of integration measure for anomalous and density-related parts 

In order to integrate over the density parts of the Gaussian like factors, we determine the change of integrate measure 
from the 'fiat' Euclidean, total self-energy S'S'g! pq (t) K to the density-related and anomalous doubled parts; we use the 
factorization of latter two different kinds of various blocks 'a = V and 'a ^ V and furthermore separate a background 
averaging functional (. . .)~(o) from the path integral (|3.92|3.93[) with the coset field degrees of freedom 

CT Vo 

d[S^ pq (t)K}=jV s ^J Yl d(5\ + (x,t)) d{5\-(x,t)) (<5A+(f,t)) 2 (<5A_(f,t)) 2 x (3.94) 

(5\+{x,t)) 2 - (<5A_(f,t)) 2 ) 2 x d(\-B D (x,t)\) smh(4\'B D {x,t)\) d(3 D {x,t)^j x 
JJ d(\c + (x,t)\) (sin(2|c + (£,t)|)) 2 dip+(x,t) x d(\c-(x,t)\) (sin(2|c_(;r,f)|)) 2 d<p-(x,t) x 



{S,t} 



x d(\e D {x,t)\) sinh(2\e D (x,t)\) d lD (x,t)\ =/V s Zk( II d(5X+(x,t)) d(5X-(x,t)) det 2 
x | d(cosh(4\'B D (x,t)\) S ) d/3 D (x,t)\ x ( J[ ^ d(cosh(2\G D (x,t)\)^ d-f D (x,t) x 



5X + (x.t) 8\-(x,t) 
5\%{x,t) 5\ 3 _(x,t) 



\ d(\c+(x,t)\ - \ sin(4|c+(a?,t)|)J dtp + (x,t) x | d(jc-(x,t)\ — \ sin(4|c_(a;,i)|)^ dip^(x,t) 



Zn[3] = ( exp 



2 R 2 H Jn 



+T 



dt 



EEO 



l-MfM tr 



p=± £ 



(3.95) 



exp { " ^ / +T ^EEt 1 " ^ 7/) ) Q tr b [ 5S *!W(*) * * E ^(*) * 

k 7/ JO p,g=± x 



l^/Wo p= ± g K 



ii n A tr 



ba 



x exp 



^ Y^y! f^maxjVx tr In 



a— 6 > 
x—x'\p=q > 
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x exp 



2 h jo 

+ f a,b,b' = l,2 



y E W(^ / .* // )(«+^ ) )^ ! i i+^(f- i ! f)(^+^ D 



(Oh 



6 a 



°V , 



(coset fields and it!? 



x exp 



x exp < i 



2frjc 



d *pE 



>*.-/ 






17 V (*^> ^P 




v 


- I £ P 



1 



4(fiV c 



dip 



^ ry p V MlmaJ^x tr In |Yj3 u + 4?) " _ (t, i' = t) 

1 _ p,g LV / x—x';p—q 



(3.96) 



2 /-+T 

max 



di dt' £ ^TV; 4df |$ U + & Vo )v\ s , Jti,,*,) U-M \ x (coset fields and a Vg 

p,q=±x,x' 



(<>) 



We hint again to following abbreviations, already specified in (|3.60I3.81[) 



(jj + a$)^ f (t p ,t' q ) = (H a p (t p ) + a$(x,t p j)6 ab S pq S s ^ S(t p - t' q ) ; 

t— l:aa' 



(3.97) 



^A,(T-\ T; f p , i'J = f; p )^ (x, t) hif{t p ,,t' qt ) f b q : b q (x', t') - h%:£,{t P , t' q ) . (3.98) 



Further simplification of above path integrals (|3.94l3.95p with coset field degrees of freedom arises from taking a separate 
saddle point approximation of the background averaging functional (|3.96|) with respect to a variation with 5<jy (x, t p ). 

This results in definite, complex- valued functions {oy (x,t p )) of space and time where the imaginary part of (cr^(x, t p )) 
has to comply with the original, already introduced, convergence generating '— ie p terms for proper convergence and 
analytic properties of Green functions. 



4 Summary and conclusion 

4.1 Transformation to the case 'I' with static disorder reduced to stationary states 



At the end of section [2] we have already specified the path integral Zi[3] (|2.31|) with the approximation to stationary states 
in frequency space which can also be attained from simplifying the two frequency dependence of the corresponding total 
disorder self-energy (a;, uj') — > 6 U>W ' S'Eg. pq (w) a t very later steps of transformations to a nonlinear sigma model 



Zi[3] =a [F[^*,t,j^J^;3] x (4.1) 

L M Mx Jo {j^) s p9=± V J V / )tF[i,*,i,\. 

Qmax = — ; < t p < +Tq ; To = T max ; (4-2) 

1 2tt 

< L0 P < +tlo ; rio = Qmax = — \ Aw = — ; K t = T max / At ; 

At J-raax 

d^> = f +n ° duj + + r° fa. = (43) 
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+n ° dco+ r+Vo+Au) duJ _ ^ ,+n p dujp 



— A. 



0+ = — alu ; fi + = il Q ; 0_ = ; 0_ = f2 + Aw ; 
T\f 2 O 

Corresponding to (|4. 1114. 41) . we can exchange the parameter fi p IT ^ p. 691) by fj.jP in p. 951) 

1 2 r )PhV 2 7 fc hv 777 

and can perform the Fourier transformation to frequency space for the anomalous doubled one-particle part according to 
following relations 

+f r+f dt r+O d0J 



dt ... = T max . . . ->■ T max - 27i - ... ; (4.6) 

Jo J- max Jo ( 7 1 max / 

J d Wp ^t p He' uAif H p {x^ p ) fe jP H= / dcdp fij^.wp) ^»e ,BAf, = (4-7) 

of, si?' 

A|VK,a; p ) = diag^.^,^); ^d u >v)) i ( 4 - 8 ) 

JK^wJ, w p ) = f - t/i CXp{ ^A t At> ~ 1 + - * e P ) J <5 gP ?7 9 _ u p) 5 s',s ; (4.9) 

^(wj.wp) = (O&VK' Wp)) T = ( - i fi cxp{ "'^ At} ' 1 + (fr r (5?) - , e p )) V ,y 9 *( W J - Wp ) ; (4.10) 

^(wj.wp) * ^. s {io q ,u v ) = (4.11) 

= dia g (A + (f -6^,u/_) ; -i^(£>_))v^',^K-^); 



pj2 ->2 

H p (x,u)p) = -hujp - i s p + — + u(x) - [1q ; Hp(x,u}p) = -hLJp—is p + — +u(x)-fj, , (4-12) 

where we additionally define doubled, shifted fields &g(wp), ^&(w' q ) with analogous boundary conditions as in relations 
(|2. 22112. 25)) for the case of dynamic disorder (cf. also relations (|2.HI2.8|) ) 

= ( £ I J) e»»**> )° ; # ^'K) = ( fz>, q (W' At \ 5 ^£)) b ■ (4-13) 

b=l 6=2 

Hence, under restriction to stationary states the total disorder self-energy 8Y,g. pq {t) simply changes to 8Yi^ pq {uj) with 
analogous changes within the the coset decomposition Sp(4)/U(2) (g> U(2) of the coset matrices T pq (x, uS) and the involved 

generators Y pq (ui), X pq (u>) so that the factorization into eigenvalues c+(x, uj), c-(x, tS) and angular parameter Cd(x,uj) 
similarly keeps its form as in (|3.25H3.27j) and (|3.36H3.40| 

S^ pq (t)K -> 5^ pq {u)k- (4.14) 
f(S,u) = exp{-r p ; b (f, W )} ; f- 1 (x , ,u)=e X p{Y p a q b (x,Gj)} ; (4.15) 



4.1 Transformation to the case T with static disorder reduced to stationary states 
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X pq (x,uj 



Y p a hx,u> 



Y, 



D;4x4 



(X,U) 



X Dvpq (x,ui 

Gd-,2x2(x,U! 



.(x, UJ 



-X*[x,w) 



\X pq {x,u) 

( \ 22 

(°) 

\ / pa 



v / ) 



-Sc d . ++ (x : oj) 5c D -+-(x,u)) 
-5cd--\ (x,oj) Scd- (x,u>) 

= P^(x,u) Y D . ixi (x,uj) P4 X4 (x,uj) ; 



- rj p X^ pq (x,uj) r] q 

-Vp X L(x,^) Vq = 



(x pq [x,w) 

( \ 22 
(°) 



6C *D; + 



(x,uj) -Sc* d . + _(x,uj) 
5c* D ., {x,uj) —5c* D . (a?, w) 



i'<i 



X D . pq (x,u>) 



-c ++ (x,u) 
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c (x, oj) 



; P 4X 4(^,w) = 



o PiU*^) 



c* ++ (x,u) 





-c*__(x,uj) 



:= \c+(x,(jj)\ e~xp{np+(x, ui)} ; c (x, u) := \c~-(x, u)\ exp{np-(x, uj)} ; 



exp 



GD;pq{x, Uj) 







-<2.D;+-(£,w) 



■ (4.16) 

4.17) 

(4.18) 

(4.19) 
(4.20) 

(4.21) 
(4.22) 



\& D {x,w)\ exp{i7 D (f,w)} 



These similarities also hold for the block diagonal density-related parts so that one can immediately state following 
relations in place of 



6Z a D a . p Jx,u>) K 



Q^' aa (x,cj) SA;,(x,u) Q a p ? q (x*,cj) 

( ( SB d . ++ (x,uj) -5B D - t+ -(x,u) 
8Bjj. + _(x,u}) -8B D . — (x,u) 



(4.23) 



\ 
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-SB d . ++ (x,uj) 5B* d . + _(x,uj) 
\ \ -5B D ; + -(x,uj) SB D . — (x,u) J pq J 

p-5X p (x.Lj) —p-5Xp (x,uj) 

5 A a (x» = diag{ '(+<5A+(£» * -<JA_ {x , ~) ) ; (-SX+(x, w) , +5X-{$, oj) ) } = SX P (x, u) K a p a p 
^ v ' „ ' J 

a— 1 a— 2 

5A p (x,cj) = diag{<5A+(x, oj) , 5X-(x,oj)} ; 5A p (S,£j) € R ; 



Q"(f,w) 



(sE%(x,oj) fj) 
-(s±f(x,oj) fj) 



(exp{i% D (x,u))y) ; Q^(x,w) = (expji S^(f,cj)| 









Q 22 / T (f,u;) = Q"(z>) 



|®_D(f,w)| exp{i^n(f, w) } ; S fl;+ _(.i,w)eC; 
Q^T^w) 5A«(f, w ) Q^ g (f, W ) ; Q a p a q (x,uj) 

Q n ' _1 (^,a;) (p- 5A p (x,w)) Q n (x» ; 
Q 22 '- 1 ^,^) (-p-<5A p (x») Q 22 (f,w); 



Q P J(x,a;) 



(4.24) 

(4.25) 
(4.26) 

(4.27) 

(4.28) 
(4.29) 
(4.30) 



However, the real self-energy variable ay ' (x, i p ) for the hermitian contact interaction takes a dependence on the difference 
of two frequencies aoj p = uj p — uj' p which therefore reduces in the limit oj p — > uj' p for stationary states to the zero frequency 
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4 SUMMARY AND CONCLUSION 



mode of a real- valued self-energy variable (x, AUlp = 0) 



A<^=0 ( o) 



(4.31) 



Consequently, we can directly convey the integration measure (|3.94[) and path integral Zu[3] (|3.95|l with background 

averaging functional (|3.96j) to the case Zj[3\ of static disorder for the restriction to stationary states; we do not outline 
the Fourier transformations of source fields and source matrices in detail and simply introduce their results for brevity 



d[5Eg m (u>)K\ =M S vk( II d(S\+(2,u>)) d(*A_(2,w)) (<5A + (f,c)) 2 (*A_(af,w) 



(4.32) 



x (((5A+(f,w)) 2 - (6\-(x,Lu)) 2 y x d(|S£,(f,w)|) sinh(4|® fl (f,a;)|) d(i D (x,uyj x 

i d(|c + (x,o;)|) (sin(2\c + (x, w)|)^ d<,o + (a:,aj) x rf(|c_(x, w)|) ^sin(2|c_(.T, w)|)^ d<fi-(x,ui) 

{x,t} 



x d(|e.o(x,<j)|) sinh(2|er>(x,o;)|) dj D (x,uj) J = .A/^g^l d(5A+(f,w)) c?(c>A_(f, w)) det 2 

■{x,t} 



6\+(x,uj) 5\ 3 _(x,u) 



\ d(cosh(4|S D (f, w )|)) dp D (2,w)\ x J|i d(cosh(2\e D (x,u;)\)) d lD { 



(X,LJ) X 



{3d} 



x ^ d(^\c + (x,io)\ ~ jSm^\c + (x,Lu)\)^J dtp + (x,L}) x i c?^|c_(x, w)| — j sin(4|c_(x,w)|)^ d(p-(x,Lj) 



Z![3] = ( exp 



jba 



TOO ~ T_ 

d ipip;x;p u ipip;x;p ™ 



x exp 



x exp 



x exp 



1 1 



" — ( ^r\ EEt 1 - tr W%M K K 



<J=± X 



a . b 



du) 



dio. 



— i V / a,o L 



(4.33) 



^^A^trm i + 6X(T-\f) (jj + ag))" 1 + (f-^/xf) + 



x—x'yp—q 



i N? f +n du du 
x exp < - ' 



o',6"=l,2 



u y imciT > y -I max 'p',q" = ± X,X 

X j T q"q' >W ) + a Vo ) s „ s , (U] q , , U q ) 



E E^^-#W^ 



«,M'=1,2 



^ T ma ,„ / r> n = + 



E E^ 

,q,q / — ± a?' 



i + MC(f - x ,f ) + ^ o } ) 1 + (/ /^;!/A / ) + £ Vo 



^(0)n-1 



-1:6a 



b"a' 



. -.(O)s -l;6"a- , „ > 

(i5 + a y J K»,^p0 



^ ^;k( w P') 



^(0) 



(coset fields and er 



^(0) 
Vb 



^(0) 

o 



(4.34) 



4.1 Transformation to the case T with static disorder reduced to stationary states 
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x exp 



t T„ 



x exp < i— 
n 



2HV -i e p J~ 

x\p—± 

K? f +{l dLU duj' 



rZj_ JC U I, T .... J ~T P>1 LV / x=x';p=q 



(r ma J (tJU p,q=±$,# 



]rr E E^ T ™*4*'K) (•^ 11 + < ) )'7 ?/ >^^)i^K) x 



(0) 



x ^coset fields and cr 

The corresponding Green function term of p. 971) and gradient operator (|3.98p has a similar form in frequency space under 
the simplifying restriction to stationary states 

+5'S ) )^(Wp,aO = (H a p ^ P ) + a^(x))s ab S pq 5 S;S , 5{lo p - u' q ) ; (4.35) 

8Kt fS ,{T-\f-u p ,u' q ) = f^'^a;)^^^,^)^^,^)-^^,^). (4.36) 

Although it is straightforward to integrate over the density-related variables d\T>r){x,oS)\, d/3u(x,LL>) and d(6A±(x, oj) J 
within the Gaussian factors of <|4.33[) in a direct manner, there occurs a further simplifying fact which concerns 

the parameter u.^ (|4.5p ; as one takes the limit ZNT t = (T max / At) — > oo for an infinite number of discrete time steps, the 

parameter fjj> approaches zero for an even simpler calculation of the density- related Gaussian factors within Zi[S] f|4.33[) . 
Considering this simplifying aspect, we finally attain the reduced path integral (|4.37j) for the case of static disorder with 
the appropriate scale of the disorder parameter and the parameter J\f x (|1.3|) for the number of spatial points which allows 
for a saddle point approximation in the limit Af x — > oo 



exp 



i_i K 

2R 



x exp { - -—J Af x 



1 1 



2 "tTT" E E tr i^^;?;p K ^inb;x;p K \ [ ^[^x^.pq ( w ) K. 
I p=± $ a ^ J 



(4.37) 



p,q=± x 



x exp 



27rf7 



p=± X 



dui 



x exp ^2lT 



2 7c. it 
% Nf duj duj 



a, 6 
a',6"=l,2 



(o)\-i 

Vb J 



(w, W = w) 
x—x'\p—q 



— e E^ T m-^S»K«)^ 



11 (tJL) (.T.JaJp'fpL± 
^(0)%-l;6'=6. 



a,M'=l,2 



o (tctt) 



E E^ 

p,q,q' = ± 



i + flKCf-Sf) (^+^ ) )- 1 + (f- 1 J f^a/^f) (s+^ o 



^(0)x-l 



-1:6a 



K,w p ) T )' aa (x,u) Uuj"„,w p ,) + 

r' f I _.. _ 



b" a 
It 



(A , — (0)\-l;6"a', // , 



^ J^-x^p') 



^(0) 
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